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In order to gain a better understanding of the Ising model in higher dimensions we have made a 
comparative study of how the boundary, open versus cyclic, of a d-dimensional simple lattice, for 
d — 1, . . . , 5, affects the behaviour of the specific heat C and its microcanonical relative, the entropy 
derivative —dS/dU. 

In dimensions 4 and 5 the boundary has a strong effect on the critical region of the model and for 
cyclic boundaries in dimension 5 we find that the model displays a quasi first order phase transition 
with a bimodal energy distribution. The latent heat decreases with increasing systems size but for 
all system sizes used in earlier papers the effect is clearly visible once a wide enough range of values 
for K is considered. 

Relations to recent rigorous results for high dimensional percolation and previous debates on 
simulation of Ising models and gauge fields are discussed. 



I. INTRODUCTION 

The Ising model on 2 and 3-dimensional lattices is 
probably the most studied model in the theory of phase 
transitions from the condensed matter point of view. 
However, since at least the early 90 's higher dimensional 
lattices have become increasingly more important due 
to their connection to gauge field theory and particle 
physics. In particular the Ising model on the 4- and 
5-dimensional cubic lattices have been studied both the- 
oretically and via simulation. In these dimensions it is 
known that in the thermodynamical limit [1, 2] the model 
follows the mean field critical exponents, but the finite 
size behaviour of the model has been hotly debated 0- 
[l5j . Part of the importance of these debates comes from 
the fact that it directly impinges on the methods used to 
derive other quantities from the finite size data, whose 
values are not already known in the limit. These meth- 
ods e.g. affect what has been done in field theoretical 
studies of the Higgs particle mass as well as other prop- 
erties studied using lattice quantum chromodynamics. 

In the work presented here we set out to make a sys- 
tematic comparison of the finite size behaviour of the 
Ising model on lattices with open and cyclic boundary 
conditions. We include both the well studied 2- and 3- 
dimensional lattices and the earlier debated lattices in di- 
mension 4 and 5. To our knowledge the previous studies 
in higher dimensions have focused exclusively on lattices 
with cyclic boundary conditions, but as recent rigorous 
results on percolation [1(1 [TtJ has shown, for this, sim- 
pler, model there are non-vanishing boundary effects in 
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higher dimension. We believe that there is a risk in fo- 
cusing only the cyclic case. 

As our sampled data reveals, there is a striking change 
in how the boundary conditions affect the finite size be- 
haviour as we pass the critical dimension D = 4. For 
D = 4 the finite size effects are more visible in the mi- 
crocanonical ensemble than in the canonical. This is 
due to the smoothing effect of fluctuations in the en- 
ergyof the system. Here the data, as previously shown 
in [18j . favours the conclusion that the specific heat at 
K c is bounded, which is consistent with the rigorous re- 
sult a — but in conflict with predictions of a weak 
divergence by renormalization theory [iljj . We also find 
that for a certain range of small sizes and open boundary 
conditions the specific heat has two maxima close to the 
critical point, one of which disappears as the system size 
increases. In the microcanonical ensemble this effect re- 
mains visible even for the largest systems we have studied 
(L=40). For cyclic boundary conditions not such effects 
are visible in either ensemble. 

For D = 5, where the specific heat is rigorously known 
to be bounded [2j], the model display more striking dif- 
ferences between the two boundary conditions. As for 
the 4-dimensional case there is a system size range with 
two maxima in the specific heat for free boundary con- 
ditions, but not for the cyclic case. This effect is again 
visible for even larger size in the microcanonical ensem- 
ble, because the smoothing effect of energy fluctuations 
is not present here. There is also a larger gap between 
the effective critical points for each size, e.g. the val- 
ues of K at which the model has the maximum specific 
heat, but this is to be expected due to the relatively large 
boundary size. However, the most striking difference is 
that for cyclic boundary conditions the model displays 
a quasi-first order phase transition. For K close to K c 
the energy distribution becomes bimodal, just as for a 
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high-g Potts model. As the system size increases the gap 
between the two maximima scales as L~ 5 / 2 , thereby giv- 
ing the model an effective latent heat for finite size, but 
not in the thermodynamic limit. Due to this bimodal be- 
haviour simulations become sensitive to small changes in 
K, as this makes the model favour one of the two maxima 
and thereby a different internal energy. In fact, some of 
the disagreements between different simulations and data 
mentioned in [3Ml5| might be due to this overlooked ef- 
fect. 



II. NOTATION 

We will investigate the behaviour of two types of lattice 
graphs, the cartesian graph product of d cycles, i.e. Cl x 
• • • x Cl, and the product of d paths Pj, X • • • X Pl for 
d = 2, 3, 4, 5. Both types of graphs have n = L d vertices, 
the product of cycles has dL d edges and the product of 
paths has dL^ 1 (L — 1) edges. They differ only at their 
boundaries where the cycle product has cyclic boundary 
conditions while the path product has open boundary 
conditions. 

We define the energy of a state a — (<7i, . . . ,a n ), with 
(7; = ±1, as E(a) — j] a i a j (summation over the 
edges {i,j}), and the magnetisation as M{a) = J2i a i 
(summation over the vertices). 

We look into two classes of quantities. The first of 
these, the combinatorial quantities from the microcanon- 
ical ensemble, provide an extremely detailed picture of 
the inner workings of the model, and they all depend 
on the energy U = E/n. The coupling K, defined 
as K(U) = —dS/dU, is of special interest here. Here 
S(U) = (1/n) In a(.E) where a(E) denotes the number of 
states a at energy E. It is described in great detail in 
[20l | how to obtain these from sampled data. 

The second class of quantities, the canonical (or physi- 
cal) ones, are the cumulants, i.e. derivatives of In Z where 
Z is the partition function. How to obtain these quanti- 
ties from the coupling function K(U) is described in (2lj . 
We will refer to them in their dimensionless forms. We 
define the internal energy as U{K) = (E) /n, the spe- 
cific heat as C(K) = var (E) jn and the susceptibility as 
x(K) = var (M) jn = ^M 2 ) jn. For each coupling K we 
receive a distribution of energies having density function 
Pr E and again we refer the reader to [2l| for how to 
obtain this distribution from the if([/)-function. 

Occasionally we will, when necessary, subscript func- 
tions with the linear order L of the lattice, as e.g. Cl(K). 
We denote by K* the coupling where Cl is at its max- 
imum. We denote by U* an energy where dK/dU has 
a local minimum and we use U + for a second minimum 
where U* <U + . 



III. 2D-LATTICES 

For 2-dimensional lattices there is no qualitative dif- 
ference between the cycle and the path product for the 
specific heat's behaviour. In Figure Q] we plot the spe- 
cific heat for both lattice types for a range of linear or- 
ders. Though the specific heat is somewhat smaller for 
the path products, not so strange considering that these 
graphs have fewer edges, the only substantial feature of 
the path products is that the peak is always located to 
the right of K c rather than to the left as it is for cycle 
products. We should here mention that for cycle prod- 
ucts with L < 320 we rely on exact data computed in (22j | 
and for L = 512 on sampled data. For path products we 
have exact data only for L < 16, see [23|], and sampled 
data otherwise. 




FIG. 1: (Colour online) Specific heat C(K) for 2D cycle prod- 
ucts (red) of linear order L = 8, 16, 32, 64, 128, 256, 512 and 
for path products (blue) of linear order L — 8, 16, 32, 64, 128, 
256. Black vertical line at K c = 0.44068 .... 

Let us compare how the maximum specific heat grows. 
We know that for square cycle products the maximum 
grows as 1.03684+ (8/ir) \nL, and we extract this from 
|24| . Note that for Cap x Cl the additive constant is 
instead 0.967550, see [25[ for an exact expression. In 
Figure [5] we show maxCL — {8/n) \nL versus 1/L for 
both types of lattices. Though we do not know what the 
limit additive constant will be for path products we can 
estimate it to be —2.30. 

The locations of the maxima are plotted in Figure [3] 
Note that the asymptotic behaviour is already known to 
be K c -K* - 0.15878/L, see (HI, for the cycle products. 
For path products there are obviously some higher order 
correction terms at work as well. 

Looking from a microcanonical perspective and plot- 
ting the derivative dK/dU, which is effectively a micro- 
canonical inverse of the specific heat, we see a more or 
less similar behaviour between the lattice types in Fig- 
ure HI This time the lattice types agree on which side of 
the asymptotic location the local minimum should be. 
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FIG. 2: (Colour online) maxC L - (8/tt) hxL vs 1/L for 2D 
cycle products (red points) of linear order L = 12, 14, 16, 20, 
24, 32, 40, 48, 56, 64, 72, 80, 96, 112, 128, 160, 256, 320, 512 
and for path products (blue points) of linear order L = 12, 
14, 16, 32, 64, 128, 256. 
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FIG. 4: (Colour online) dK/dU for 2D cycle products (red) 
of linear order L = 16, 32, 64, 128, 256, 512 and for path 
products (blue) of linear order L — 16, 32, 64, 128, 256. Black 
line at U c = v2. 
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FIG. 3: (Colour online) Location K* of maximum specific 
heat vs 1/L for 2D cycle products (red points) of linear order 
L = 16, 20, 24, 32, 40, 48, 56, 64, 72, 80, 96, 112, 128, 
160, 256, 320, 512 and for path products (blue points) of 
linear order L = 16, 32, 64, 128, 256. The black asymptote is 
K c - 0.15878/L. 



FIG. 5: (Colour online) min<9L" / 'dll for 2D cycle products 
(red points) of linear order L = 12, 14, 16, 20, 24, 32, 40, 48, 
56, 64, 72, 80, 96, 112, 128, 160, 256, 320, 512 and for path 
products (blue points) of linear order L = 12, 14, 16, 32, 64, 
128, 256. Black line is the asymptote vr/(8 InL). 



For these lattices we can actually give the asymptot- 
ical behaviour for how the minimum of dK/dU should 
approach 0. Since the maximum specific heat grows 
as {8/tt) InL, and we excpect no asymyptotic difference 
here between the lattice types, the minimum dK/dU 
must vanish at exactly the inverse rate 7r/(8 InL). This 
is shown in Figure [3] Note that the path product lat- 
tices approach the asymptote from above and the cycle 
products from below. 

IV. 3D-LATTICES 

For 3-dimensional lattices we see no major difference 
in behaviour between the two lattices types either. In 



three dimensions we do not an exact solution to rely 
on, so we can not discern as many details as for 2D- 
lattices. For cycle products we have sampled data for 
L = 6, 8, 12, 16, 32, 64, 128, 256, 512, as used in [26] while 
for path products we have sampled data only for smaller 
lattices; L = 6,8,12,16,24,32,46,64,96,128. For both 
lattice types we use exact data for L = 4. 

Again let us compare the lattice types. The specific 
heat is shown in Figure [6j As in the 2D-case the path 
products' maximum specific heat is smaller. The location 
of the maximum approach K c from above for both lattice 
types, see Figure [7] However, for path products this 
maximum is located much farther away from K c . The 
black curves in the plot were obtained from a best fit of 
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all available data points to the simple scaling formula 

K* = c Q + a L- Xl + c 2 L- X2 (1) 

where we pre-ordained c$ = 0.2216546 as found in [26[. 
For the cycle products we received c\ = 0.3775, Ai = 
1.654, c 2 = -0.6169 and A 2 = 2.297. For the path prod- 
ucts this gave c\ = 2.292, Ai = 1.617, c 2 = -2.101 and 
A 2 = 1.961. The fitted curves agree beautifully with the 
points for all L and they are included in the plot. 



values for both lattice types. Note that for cycle products 
this minimum is actually increasing for L < 16, whereas 
for path products they are consistently decreasing. This 
decrement speeds up dramatically, as it should, for L > 
64 as shown in the figure. 
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FIG. 8: (Colour online) dK/dU for 3D cycle products (red) 
of linear order L = 6, 8, 12, 16, 32, 64, 128, 256, 512 and for 
path products (blue) of linear order L = 6, 8, 12, 16, 24, 32, 
48, 64, 96, 128. Black line at U c w 0.99063. 



FIG. 6: (Colour online) Specific heat C(K) for 3D cycle prod- 
ucts (red) of linear order L = 4, 6, 8, 12, 16, 32, 64, 128, 
256, 512 and for path products (blue) of linear order L = 4, 
6, 8, 12, 16, 24, 32, 48, 64, 96, 128. Black vertical line at 
K c w 0.2216546. 
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FIG. 9: (Colour online) Minimum dK/dU for 3D cycle prod- 
ucts (red) of linear order L = 6, 8, 12, 16, 32, 64, 128, 256, 
512 and for path products (blue) of linear order L = 6, 8, 12, 
16, 24, 32, 48, 64, 96, 128. 



FIG. 7: (Colour online) Location K* of maximum specific 
heat for 3D cycle products (red) of linear order L = 4, 6, 8, 
12, 16, 32, 64, 128, 256, 512 and for path products (blue) of 
linear order L = 4, 6, 8, 12, 16, 24, 32, 48, 64, 96, 128. See 
text for parameters of fitted curves. 

On the microcanonical side we note that the derivative 
dK/dU is much smaller for cycle products than for path 
products, see Figure[5] In Figure|H]we show the minimum 



V. 4D-LATTICES 

In [l8| the authors studied the Ising model on cycle 
products and, unexpectedly, found that the specific heat 
seemed to be bounded at K c , a result which is compatible 
with the rigorous bounds but in disagreement with non- 
rigorous renormalization predictions. As we shall now 
see, unlike in the lower dimensional cases, in 4 dimen- 
sions there are also distinct differences in the finite size 
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behaviour of the cycle and path products. In short, the 
specific heat of the path and cycle products differ some- 
what in their behaviour, however on the microcanonical 
side larger differences are clearly visible. The minimum 
of dK/dU at U* is much higher for path products and is 
decreasing with L, compared to the cycle products where 
it is actually increasing. Also, for L > 12 there is a sec- 
ond local minimum in dK/dU at U + , where U* < U + . 
This minimum becomes a global minimum for L > 24. 
We will now proceed as for the other lattices. 

The specific heat is shown in Figure [TO] Note espe- 
cially the peculiar behaviour of C near K* for the path 
products. This is a phenomenon due to the extra local 
minimum at U + in dK/dU . Due to this peculiarity the 
location of the global maximum will display an irregular 
behaviour, depicted in Figure [TT1 



o.no - 



0.165 - 



0.160 



0.155 - 



0.150 




do 64 32 24 16 12 



1/L 



FIG. 11: (Colour online) Location K* of the maximum spe- 
cific heat C(K) for 4D cycle products (red) of linear order 
L = 4, 6, 8, 10, 12, 16, 20, 24, 32, 40, 48, 56, 60, 64, 80 and 
for path products (blue) of linear order L = 4, 6, 8, 12, 16, 
24, 32, 40. 
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FIG. 10: (Colour online) Specific heat C(K) for 4D cycle prod- 
ucts (red) of linear order L = 4, 6, 8, 12, 16, 24, 32, 40, 48, 
64, 80 and for path products (blue) of linear order L = 4, 6, 
8, 12, 16, 24, 32, 40. Black vertical line at K c w 0.1496947. 
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In Figure[TT]we have fitted curves to K* that are forced 
to have the same Co = K c = 0.1496947, obtained in [l8j . 
and obviously there is no conflict in this. However, the 
correct asymptotic behaviour for the path products is 
only obtained from L > 16, for our choice of lattice sizes, 
leaving us only four points to fit our curve to. Thus we set 
C2 = in both cases and obtain C\ = 0.243 and Ai = 2.03 
for the cycle products. For the path products we obtain 
c\ = 0.166 and Ai = 0.962. Strikingly, the exponents 
differ by about a factor two. 

Let us now turn to the microcanonical side. In Fig- 
ure [12] we show dK/dU for both lattice types, near U c , 
where each curve has a local minimum. Outside the 
picture there also materialises, for path products with 
L > 12, a second local minimum, see Figure 1131 This 
minimum is not as sharp as the other minimum, thus 
making it difficult to pin-point it to a high accuracy. 
Moreover, this minimum becomes a global minimum for 
L > 24. This phenomenon is shown in Figure IT4l 

Here the picture gets more complicated. First, as was 
demonstrated in [18|], the values of the local (and only) 



FIG. 12: (Colour online) dK/dU for 4D cycle products (red) 
of linear order L = 4, 6, 8, 12, 16, 24, 32, 40, 48, 64, 80 and 
for path products (blue) of linear order L = 4, 6, 8, 12, 16, 
24, 32, 40. Black line at U c « 0.770527. 



minima for the cycle products are increasing to about 
0.0042, thus leading to an asymptotically finite specific 
heat. As discussed in discussed in [l8[ this is in dis- 
agreement with a prediction from renormalization theory, 
which predicts a logarithmic singularity in the specific 
heat for D — 4. It is rigorously known Q that the cirit- 
cal exponent a — for D = 4. In [l8| data from both a 
single spin Metropolis and a Wolff-cluster algorithm [27j 
were compared to rule out methodological errors and the 
conclusion was that either systems as large as L = 80 are 
still dominated by finite size effects or the renormaliza- 
tion prediction fails for the D = 4 case, which is known to 
be the critical dimension for the nearest neighbour Ising 
model. 

Second, the local minimum nearest U c for the path 
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FIG. 13: (Colour online) dK/dU for 4D path products of 
linear order L = 4, 6, 8, 12, 16, 24, 32, 40. 
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FIG. 14: (Colour online) Minimum value of dK/dU for 4D 
cycle products (red, converging to 0.004) of linear order L = 4, 
6, 8, 10, 12, 16, 20, 24, 32, 40, 48, 56, 60, 64, 80 and for path 
products (blue) of linear order L = 4, 6, 8, 12, 16, 24, 32, 40. 
The blue points converging to 0.019 are the minima located 
at U* and those converging to 0.009 are located at U + . 



products seems to approach a different and even greater 
vafue of roughly 0.019. Third, the other local minimum 
for the path products, also becoming a global minimum, 
could approach yet another value of about 0.009 in be- 
tween the other two. Forcing this minimum to approach 
the same minimum as for the cycle products seems out 
of the question with our current set of data though. 

In Figure [14] the fitted curves a la ([1]) use c-i = and we 
received cq — 0.0042, c% = —0.016, when setting Ai = 1 
for the cycle products, see a discussion on this in [l8| . 
For the first set of local minima for the path products, 
the ones nearest U c we received cq = 0.0190, ci = 0.223, 
Ai = 2.47. The second set for the path products, farther 
away from U c , gave us Co = 0.00870, ci = 0.0568, Ai = 
0.583 when the curve was fitted to all data points. 

In Figure [15] we show the location of the local minima 



FIG. 15: (Colour online) Location U* of the minimum value 
of dK/dU for 4D cycle products (red points) of linear order 
L = 4, 6, 8, 10, 12, 16, 20, 24, 32, 40, 48, 56, 60, 64, 80 and for 
path products (blue points) locations U* and U + for linear 
order L = 4, 6, 8, 12, 16, 24, 32, 40. 



for both graph types. The location of the minima for the 
cycle products are clearly approaching U c m 0.770527, 
see [13|. However, the local minimum U* for the path 
products does not appear to approach the same limit for 
our current set of points. We believe this is only due 
to the lattices being too small; our current record size 
is only L = 40 for the path products. The second min- 
imum's location U + could very well be the same as for 
the cycle products limit minimum location. An alterna- 
tive possible scenario is that the local minimum U* will 
indeed have a different limit than U c , but that this does 
not effect the behaviour of the canonical quantities since 
their behaviour is guided by the global minimum, the 
location of which appears to approach U c . 

The location U* of the minimum for the cycle products 
were set to have the limit cq — U c — 0.770527. This 
gave ci = 4.14 and Ai = 1.77 when setting C2 = and 
fitting to the ten largest lattice sizes. This looks very 
convincing. For the path products' local minimum U* we 
received cq = 0.7528, i.e. a different limit, and c\ = 23.1, 
Ai = 2.39 when using C2 = and fitting to the last five 
data points. The local minimum U + , which is a global 
minimum for L > 24, gave an acceptably good fit using 
c = U c , ci = 0.788, Ai = 0.355. 



VI. 5D-LATTICES 

For the 5-dimensional lattice the difference in be- 
haviour between path products and cycle products be- 
comes even more pronounced. Let us begin with the spe- 
cific heat as usual. As Figure [TBI shows . the behaviour is 
similar to that of the 4-dimcnsional lattice. 

For these lattices the specific heat is rigorously known 
[2J to be bounded. Using our scaling formula UJ) with 
c 2 = we received cq = 172.8, ci = —383.1 and 
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FIG. 16: (Colour online) Specific heat C(K) for 5D cycle prod- 
ucts (red) of linear order L = 4, 6, 8, 10, 12, 16, 20, 24 and 
for path products (blue) of linear order L = 4, 6, 8, 12, 16, 
20. Black vertical line at K c w 0.113914. 



FIG. 18: (Colour online) Location K* of the maximum spe- 
cific heat C(K) for 5D cycle products (red) of linear order 
L = 4, 6, 8, 10, 12, 16, 20, 24 and for path products (blue) of 
linear order L = 4, 6, 8, 12, 16, 20. 



Ai = 1.434 giving an excellent fit over all the lattices. 
This is demonstrated by the red points and their fitted 
curve in Figure [T7] Ideally we should also receive the 
same asymptotic value for the path products though the 
finite size scaling may be quite different. In the figure we 
have forced cq to the same value, 172.8, as for the cycle 
products. This gave a = -173.1 and Ai = 0.2222. The 
fit does not look too strained, but forcing the value of cq 
was necessary for this. 

max C 
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FIG. 17: (Colour online) Maximum specific heat ma,xC(K) 
for 5D cycle products (red) of linear order L — 4, 6, 8, 10, 12, 
16, 20, 24 and for path products (blue) of linear order L — 4, 
6, 8, 12, 16, 20. 

The locations of the maximum specific heat behave in 
the typical smooth fashion for the cycle products but the 
path products have the peculiar behaviour found in the 
4D-lattice. Here, for L > 8 the maximum has jumped 
farther out from K c . This shows clearly in Figure [18] 
where we plot the location of the maxima. 



For the fitted curves we have used cq — K c = 0.113914, 
see [l]|. Our best fit gave ci = 0.183 and Ai = 2.44 for 
the cycle products, but asymptotically we should have 
Ai = 5/2. For the path products we found c\ = 0.135 
and Ai = 0.918, based on the points for L > 8. The 
exponents thus differ considerably so there could be a 
different scaling rule in play here. 

Moving on to the microcanonical side we show in Fig- 
ure[19]the dK/dU-cmves for the two types of lattices. As 
for the 4D-lattices we again see a huge qualitative differ- 
ence between them. To begin with, for the path products 
we see a second local minimum in the curve, located far 
away from U c , which becomes the global minimum for 
L > 12, see Figure [501 Also, the minima of the cycle 
products are all negative, whereas for the path products 
they are all positive. This is a striking difference, since 
a negative minimum is normally seen as an indication of 
a first-order phase transition, and as we shall see later 
there is indeed a first-order-like behaviour at the critical 
point for the cycle products. 

As for the values of the minima we turn our eyes to 
Figure [2TJ The fitted curve for the cycle products has 
been forced to Co = 0, which may or may not be correct. 
Then ci took the best-fit value —0.0054 and Ai chose 
0.447, when setting also c 2 = 0. The point for L = 20 
was left out of the fitting process since it clearly deviates 
from the others in its behaviour. 

For the path products we now have two sets of minima 
to look at, just as for the 4D-case. We run into trouble 
for the minima nearest U c though. Fitting (TTJ) with only 
one exponent gave negative Ai. This would imply that 
the minimum goes to infinity, which in its turn implies 
that the specific heat would go to zero. This is clearly 
an inconsistent behaviour due a too simple fit. Using 
two exponents instead, a best-fit gave cq = 0.062, c\ = 
-0.071, c 2 = 0.042, Ai = 0.151 and A 2 = 1.12. The 
fit is very good. However, it should be kept in mind 
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FIG. 19: (Colour online) dK/dU for 5D cycle products (red) 
of linear order L = 4, 6, 8, 10, 12, 16, 20, 24 and for path 
products (blue) of linear order L — 4, 6, 8, 12, 16, 20. Black 
line at U c « 0.67549. 



0.018 - 




0.016 - 



0.014 - 



0.012 - 



0.010 



0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 



FIG. 20: (Colour online) dK/dU for 5D path products of 



that this means we have fitted 5 parameters to 6 points, 
so the value of Co (as the other parameters) is likely to 
change upon receiving further data for larger lattices. 
What seems clear is that this minima has an increasing 
trend with increasing L. 

Considering instead the minima farther out from U c we 
see in the figure that it becomes a global minimum for 
L > 12. Also, this local minimum only exists for L > 8. 
The fitted curve, setting C2 = 0, has the parameters cq 
0.007.'). ci = 0.044, Ai = 0.77, which gives a convincing 
fit but it relies on only 4 points. 

Figure displays the location U* of the minimum 
dK/dU- values. For the path products we thus see two 
sets of these, U* and U + . The curve fitted to the 
red points, i.e. the cycle products, has been forced to 
cq = U c = 0.67549 and this gave a = 5.67 and A x = 2.34. 
Setting Co to U c for the path products' minima at U* gave 
the curve a very strained look. A free fit instead gave 



FIG. 21: (Colour online) Minimum value of dK/dU for 5D 
cycle products (red points) of linear order L = 4, 6, 8, 10, 12, 
16, 20, 24 and for path products (two sets of blue points) of 
linear order L = 4, 6, 8, 12, 16, 20. The blue points converging 
to 0.06 (outside the plot) are the minima at U* and the blue 
points converging to 0.007 are the minima at U + . 



c = 0.6478, ci = 13.6 and Ai = 2.67. Receiving dif- 
ferent limit values (which we also did for the 4D-lattice) 
is a rather unsatisfactory outcome, but the data points 
are very nicely fitted to this curve, and looks good even 
if we were to zoom in on these points only. More likely, 
the lattices are far too small to have received a correct 
asymptotic behaviour. However, we should not rule out 
the possibility of a second asymptotic critical point. Af- 
ter all, what matters to the model is that the global min- 
imum goes to the same value as for the cycle products. 
This minimum's location U + is also seen in the figure but 
it is very far from U c . The curve through the last three 
points only, has a forced c = U c = 0.67549 and this gave 
d = 0.946 and Ai = 0.218. 

One focus of the earlier mentioned debates regarding 
the 5-dimensional Ising model (IH [l3[ was the scaling 
with L of the susceptibility at K c . Here the main term in 
the scaling has been predicted (f| [28[ to be L 5 / 2 for cyclic 
boundary conditions, whereas an ordinary scaling would 
have been proportional to L 1 / 



L 2 . In addition (11 1 
predicted that there should be additional terms leading 
to a maximum in at L = 9 after which it should 

decrease to a non-zero limit, a prediction which was crit- 
icized in [lj|. Given our current, larger, data set we have 
made a brief return to this question, for both boundary 
conditions. In Figure[22]we see xjl^l 2 for several values 
of K close to K c . That the quotient is converging to a 
finite non-zero limit seems quite plausible for a range of 
possible K c , but with the limited system sizes available 
to us it is hard to say what the limit should be. In [ll], [l3[ 
data for K = 0.1139150 were used, corresponding to the 
orange points in our figure. However, this value is larger 
than our best estimate and as the figure shows the quo- 
tient is extremely sensitive to even small changes in the 
value of K c . The data does support a maximum in the 
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quotient, but a location at L = 12 seems more likely than 
the L — 9 predicted in [11]. 

In Figure[231 we plot the quotient xjL 2 , corresponding 
to x/L 1/v , for open boundary conditions. We see an 
excellent agreement over the full range of system sizes. 
We thus conclude that for open boundary conditions the 
dominant term in the scaling is L 2 and that any cor- 
rection terms are significantly smaller. Experiments also 
showed that this quotient was much less sensitive to the 
value of K c than for cyclic boundaries. 

X/L 5 ' 2 
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FIG. 22: (Colour online)x/£ 5/2 vs L for 5D cycle products 
of linear order L = 4, 6, 8, 10, 12, 16, 20, 24 at K = 0.1139139 
(red points), K = 0.1139150 (orange points), K = 0.1139155 
(green points). The purple lines correspond to the upper and 
lower ends of the estimate K c = 0.1139139(5). 
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FIG. 23: (Colour online)x/£ 2 vs 1/L for 5D path products 
of linear order L = 4,6, 8, 12, 16, 20 at K = 0.1139139 (blue 
points). The line has slope 0.825 and the asymptotic value is 
0.0821. 




24 20 16 



1/L 



FIG. 24: (Colour online) Location U* of the minimum value 
of dK/dU for 5D cycle products (red points) of linear order 
L = 4, 6, 8, 10, 12, 16, 20, 24 and for path products (U* and 
U + ) of linear order L = 4, 6, 8, 12, 16, 20. The points at U + 
are above those at U* . 



A. The nature of the phase transition for the cycle 
products 

The 5-dimensional cycle products are quite exceptional 
since the dK/dU become negative. This is then obviously 
reflected in the function K(U) by having a local maxi- 
mum and a local minimum. This, in turn, means that the 
distribution of energies becomes bimodal for a range of 
temperatures, see [2l|, [29[ for discussions from the math- 
ematical and computational perspectives respectively. A 
bimodal energy distribution is normally the clearest sign 
of first-order, or discontinuous, phase transition, with a 
latent heat given by the distance between the two max- 
ima of the energy distrubution. For D > 4 the Ising 
model has been rigorously proven [l|,|3(| to have a second- 
order phase transition with critical exponents given by 
the mean field model. Hence we here see a clear differ- 
ence between the asymptotic and finite size behaviour, 
which we will now examine closer. 

In Figure [23] we show the normalised distributions at 
the temperature where the two distribution peaks are 
equal and the inset picture shows a zoomed-in plot clearly 
demonstrating that the distributions are indeed bimodal. 
Denoting the location of the local maxima of this bimodal 
distribution as U\ and U2 their difference goes to zero 
as roughly U2 — U\ ~ 6.2/L 5 / 2 . The right maximum U2 
behaves as roughly U c + 13/L 5 / 2 . The temperature where 
the maxima of the distribution are of equal height scales 
as K c + 0.12/L 5 / 2 . 

These results mean that for finite L we see an effective 
latent heat with size proportional to L -5 / 2 and, as for 
other models with a first order phase transition in their 
thermodynamic limit, the finite size model demonstrates 
a bistable behaviour at K c . This behaviour is noticeable 
in the simulations, where the system will spend a long 
time at one maximum and then rapidly switch to the 
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other. 

Given the scaling of XJ% — U\ we find it quite possible 
that the model will display this bimodal distribution for 
all finite L and only reach the pure mean field behaviour 
in the thermodynamic limit L — » oo. 
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FIG. 25: (Colour online) Normalised distribution of energies 
at the temperature where the two peaks of the density func- 
tion are equal, for 5D cycle products of linear order L = 4, 
6, 8, 10, 12, 16, 20, 24 (blue). The inset picture shows a 
zoomed-in version. 



rather non-spectacular behaviour near K c as Figure [27J 
shows. Even at the temperature where the specific heat 
is at its maximum the distribution stays quite gaussian. 
However, at the temperature range where specific heat 
grows most quickly, there is some typical behaviour. For 
example, the usual left-to-right skewness is found, but it 
is considerably less pronounced than for the cycle prod- 
ucts. 
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FIG. 27: (Colour online) Density functions for the 5D path 
products L = 16 near K c at K = 0.1164, 0.1166, 0.1168, 
0.1170, 0.1172 (left to right). 



The density function of the energy distribution for the 
cycle products undergoes some dramatic shape-shifting 
near K c . It starts as a gaussian function, then becomes 
heavily skewed to the left, then bimodal, then skewed 
to the right, and finally it becomes gaussian again. In 
Figure [26] we show this as a gallery of distributions for 
L = 16 near K c . 
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0.68 0.70 0.72 0.74 

FIG. 26: (Colour online) Density functions for the 5D cycle 
products L — 16 near K c at K = 0.11396, 0.11398, 0.11401, 
0.114027, 0.11405, 0.11410, 0.11415 (left to right). 

The path products on the other hand, for which 
dK/dU stays positive, clearly can not have a bimodal 
distribution. Consequently, the density function has a 



VII. CONCLUSIONS 

As we have seen, the finite size behaviour of the Ising 
model becomes significantly more complex as we pass 
the critical dimension D = 4. In the case of 4 dimensions 
we have found that an additional maximum in the spe- 
cific heat appears for a certain range of sizes, but only 
for open boundary conditions. We also found that both 
the cyclic and open cases seem to favour a bounded spe- 
cific heat, with a value which they approach from below 
and above respectively, thereby contradicting a predic- 
tion from renormalization theory [l9[ . 

Some of the peculiarities of the 4-dimensional cases 
appears in the 5-dimensional model as well, but here the 
most striking feature is the quasi first order phase transi- 
tion for cyclic boundary conditions. The presence of this, 
for the Ising model, unusual phase transition type, has 
numerous implications for Monte Carlo experiment and 
their analysis. 

1. The model is highly sensitive to small changes in K. 
In Figure [2H the two curves closest to the bimodal 
one corresponds to the values K = 0.11401 and 
K = 0.11405. Here changes in the 5th decimal 
gives a significant change in the expected internal 
energy of the model. 

2. Simulations at values K at the balanced bimodal 
value will display a meta-stable behaviour, in which 
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the models will stay for some time at one maximum 
and then rapidly change to the other. 

3. Unlike what was initially hoped, the Wolff-cluster 
algorithm 27] suffers from critical slowing down 
at first-order phase transitions, as was rigorously 
proven in [31(. This means that, just as for 
the Metropolis algorithm, additional care must be 
taken in the study of this model in order to make 
sure that the generated data does not suffer from 
unwanted correlations 

4. If only the moments of the energy distribution are 
studied at a given temperature the clearly non- 
gaussian form of the energy distribution can poten- 
tially lead to incorrect conclusions regarding higher 
moments based on the low order moments. There 
will e.g. be significantly non-zero odd moments for 
a longer interval of K near K c . 

As we mentioned in the introduction the finite size be- 
haviour of the simple percolation model has also been 
studied in higher dimensions. In [16j Aizenman found 
that the size of the largest cluster at p c scales, at most, 
as L 4 and conjectured that for cyclic boundary condition 
the scaling should be of order L 2d / 3 instead. Recently 
this conjecture was proven in [l?} , apart from a logarith- 
mic factor which was removed in a follow up paper by 
the same authors (32j. 

The general assumption in simulation physics is that as 
the system grows larger the effect of the boundary should 
diminish and the model converge to the same thermo- 
dynamic limit independently of the choice of boundary 
condition and the rough shape of the finite region. Fol- 
lowing the seminal work of van Hove [33] we know that 
for a large class of models this will be the case, but it is 
worth to notice that van Hove's proof does in fact not ap- 
ply to a sequence of larger and larger lattices with cyclic 
boundary conditions, as the proof assumes that the edge 
and vertices of the smaller lattices are subsets of those of 
the larger. However, for K well away from K c we see no 
reason to expect anything else. Nonetheless, as the per- 
colation results indicate, it is quite plausible that when 



we stay at, or near, K c we will see drastically different 
scalings depending on the choice of boundary conditions. 

Generalizing the conjecture of Aizenman [l6| we pose 
the following conjecture for the q > 1 random cluster 
model [H 

Conjecture 1 There exists a constant c(q) and a di- 
menison D(q) such that for the random cluster model 
with q > 1 and d > D(q), the largest cluster in the 
d-dimensional cycle product with side L contains L c ^ d 
vertices at p c {q)- 

The case 9=1 corresponds to percolation and q = 2 to 
the Ising model. We intend to return to this conjecture 
in an upcoming paper. 

For models more complicated than the Ising model, 
such as those from QCD 0, M, EH, quantitative differ- 
ences between different boundary conditions have cer- 
tainly been expected. The unexpected appearance of 
qualitative differences in the nature of the finite-size 
phase transition, such as what we have seen for D = 5, 
poses a more serious problem. Close to a phase transition 
point different boundary conditions can potentially lead 
to quite distinct asymptotics, thereby making it neces- 
sary to study more than one kind of boundary in order 
to help pinpoint the correct infinite system asymptotics. 

On a more general note we believe that the sensitivity 
to K displayed in Figure [2H clearly demonstrates why 
single value simulations of models no longer suffice for 
studying models at the level of accuracy needed in mod- 
ern physics. Instead methods allowing a continuous vari- 
ation of the system parameters after the simulation, such 
that of [20| . are needed. The authors of [l3[ stated that 
"for a full understanding of the problem a variation of 
parameters over a broad range is desirable" and looking 
at our current findings we most emphatically agree. 
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